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In 2011, a discrepancy between the values of the Planck constant measured by counting Si atoms and by 
comparing mechanical and electrical powers prompted a review, among others, of the measurement of the 
spacing of ^®Si {220} lattice planes, either to conhrm the measured value and its uncertainty or to identify 
errors. This exercise confirmed the result of the previous measurement and yields the additional value 
^220 = 192014711.98(34) am having a reduced uncertainty. 
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I. INTRODUCTION 

efforts are in progress on accurate determinations of 
the Planck, h, and Avogadro, TVa, constants>i They are 
prompted by the proposal of a new kilogram definition 
based on a conventional value of the Planck constant, Na 
and h being linked by the molar Planck constant, Aa/i, 
which can be accurately measured)^ 

The most accurate way to determine Na is by count¬ 
ing the atoms in a single-crystal Si ball highly enriched 
with A measurement was completed in 2011; the 

uncertainty associated to the determination of h from 
this measurement is 3.0 x 10“®/i. The measured value 
differed from the result of a watt-balance comparison of 
mechanical and electrical powers by 3.7 times the com¬ 
bined standard uncertainty.— 

Although subsequent watt-balance determinations^ 
gave h values in substantial agreement with that obtained 
by atom counting, this prompted a reassessment of the 
Na uncertainty to identify whether errors were done. All 
the necessary measurements are being scrutinized and re¬ 
peated aiming at a smaller uncertainty, thus carrying out 
stress tests of all the technologies to confirm that the in¬ 
tended performances are met. In the present paper we 
report about the measurement of the lattice parameter by 
means of combined x-ray and optical interferometry and 
give an additional result having a reduced uncertainty. 

II. Na measurement 


where m and V are the crystal mass and volume, M is 
the mean molar mass, a®/8 is the atom volume, and 8 
is the number of atoms in the cubic unit cell. Since the 
binding energy of the Si atoms is about 5 eV and the 
mass of a Si atom is about 26 GeV, M and m can be 
viewed as the molar mass and mass of an ensemble of 
free atoms. To make the kilogram redefinition possible, 
the targeted accuracy of the measurement is 2 x 10 “®Aa. 

From ([T|) , it follows that the Na determination requires 
the measurement of i) the lattice parameter - by com¬ 
bined x-ray and optical interferometry^, ii) the amount 
of substance fraction of the Si isotopes and, then, of the 
molar mass - by absolute mass-spectrometry^—, and iii) 
the mass and volume of nearly perfect crystal-ball having 
about 93 mm diameter— rii 

Silicon crystals may contain chemical impurities, in¬ 
terstitial atoms, and vacancies, which implies that the 
measured mass value does not correspond to that of an 
ideal Si crystal and that the crystal lattice may be dis¬ 
torted. This means that crystals must be characterized 
both structurally and chemically, so that the appropri¬ 
ate corrections are applied— The mass, thickness and 
chemical composition of the oxide layer covering the ball 
must be taken into account; they are measured by optical 
and x-ray spectroscopy and reflectometry— 


III. LATTICE PARAMETER MEASUREMENT 
A. X-ray/optical interferometry 


The Na value is obtained from measurements of the 
molar volume, VM/m, and lattice parameter, a, of a 
perfect and chemically pure silicon single-crystal. In a 
formula. 


Na 


8MV 

a'^m 


( 1 ) 
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The combined x-ray and optical interferometer used 
to measure the lattice parameter is described by Fer- 
roglio et. al— As shown in Fig. [1] it consists of three 
blades, 1.20 mm thick, so cut that the {220} planes are 
orthogonal to the blade surfaces. X rays from a (10 x 1) 
mm^ Mo Kai line source are split by the first crystal 
and recombined, via a transmission crystal, by the third, 
called analyser. The interference pattern is imaged onto 
a multianode photomultiplier tube through a pile of eight 
Nal(Tl) scintillator crystals. The photomultiplier image 
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FIG. 1. Combined x-ray and optical interferometer. The yel¬ 
low line indicates the continuation of the laser beam, at 21 
mm from the analyser base, where the spacing of the diffract¬ 
ing planes was surveyed. The optical interferometer and fixed 
Si crystal rest on a common silicon plate (not shown). The 
analyser displacement and attitnde (pitch and yaw angles) 
are optically sensed via quadrant detection of the interfer¬ 
ence pattern. The transverse, y and z, displacements and roll 
rotation of the analyser are sensed via a reference 90° trihe¬ 
dron (resting on the same platform as the analyser) and three 
capacitive sensors faced to it (not shown). 


projected on the analyser is (1 x 11.2) mm^, with a pixel 
size of (1 X 1.4) mm^. 

When the analyser is moved along a direction orthogo¬ 
nal to the {220} planes, a periodic variation in the trans¬ 
mitted and diffracted x-ray intensities is observed, the 
period being the diffracting-plane spacing. The move¬ 
ment, up to 5 cm, requires to control the analyser atti¬ 
tude to within nanoradians and vibrations and position 
to within picometers. The analyser displacement and ro¬ 
tation are measured by optical interferometry; the neces¬ 
sary picometre and nanoradian resolutions are achieved 
by phase modulation, polarization encoding, and quad¬ 
rant detection of the fringe phase. To eliminate the ad¬ 
verse influence of the refractive index of air and to ensure 
millikelvin temperature uniformity and stability, the ap¬ 
paratus is hosted in a thermo-vacuum chamber. 

We measured the lattice parameter of a number of nat¬ 
ural Si crystals; the link between the results of these mea¬ 
surements and the measured value for the enriched crys¬ 
tal used to determine TVa is given by Massa et. 

All the past measurements relied on the same optical 
interferometer, that served us since 1994. In order to 
exclude systematic effects, we assembled a new one and 
integrated it in the apparatus. The main novelties of the 
upgraded system are listed hereinbelow. 

A 532 nm frequency-doubled Nd:YAG laser substi¬ 
tuted for the previous 633 nm diode laser ~ stabilised 
by frequency-offset technique against the frequency of an 
He-Ne laser which, in turn, was stabilized against com¬ 
ponent flis of the transition 11-5 R(127). The laser 
was better collimated, thus halving the correction for 
diffraction effects. Furthermore, the residual pressure in 


the vacuum chamber has been reduced by an order of 
magnitude, below 0.01 Pa. This makes any correction 
for the refractive index of the residual gas in the vac¬ 
uum chamber inessential and ensures the calibration of 
the optical interferometer with a negligible uncertainty. 

Contrary to our past measurement, the lattice spac¬ 
ing was surveyed along an horizontal line at 21 mm from 
the analyser base (instead of the previous 26 mm) and 
the correction for the self-weigh deformation was recal¬ 
culated. 

A new optical bench is clamped to the vacuum cham¬ 
ber; it collimates the laser beam, modulates the phase 
of the TT-polarized component, and delivers it to the in¬ 
terferometer by a pointing mirror and a window of the 
vacuum chamber. The delivery, collimation, modulation, 
and pointing systems - optical fiber, beam collimator and 
polarizer, phase modulator, and injection mirror - have 
been rebuild to conform to the new wavelength. 

Previously, the orthogonality between the laser beam 
and the analyser was only occasionally checked. This was 
done by observing simultaneously, via a visual autocol¬ 
limator placed - when necessary - outside the vacuum 
chamber, the analyser and laser beam through the out¬ 
put port of the interferometer. To gain the on-line control 
of the beam pointing, an home-made telescope picks up 
part of the beam delivered to the detector. In order to 
ensure stability, it is clamped on the same base plate as 
the x-ray/optical interferometer. 

A plate beam-splitter was manufactured ad-hoc and 
substitutes for the cube beam-splitter previously used 
to ensure that the difference of the transmitted- and 
reflected-light paths is insensitive to the beam transla¬ 
tions and rotations. Therefore, the components of the op¬ 
tical interferometer - beam splitter, quarter-wave plates, 
and fixed mirror - were replaced and assembled anew. 

In order to make the interfering beams parallel, the 
components of the optical interferometer are cemented on 
a glass plate supported by three piezoelectric actuators. 
As shown in Fig. [21 a noise at a frequency of about 
1 mHz caused phase instabilities between the x-ray and 
optical fringes. A new power supply was realized, having 
a sub part-per-million stability over the time scales, from 
1 s to 1 h, relevant to the lattice parameter measurement. 
Eventually, the phase noise was reduced to the shot noise 
limit of the x-ray photon count. 

The Physikalische Technische Bundesanstalt found a 
contamination of the surfaces of the x-ray interferometer 
by Cu, Fe, Zn, Pb, and Ca caused by the wet etching 
used by the INRIM to remove any residual stress due to 
surface damage after the crystal machining. The contam¬ 
ination was removed by cleaning the crystal in aqueous 
solutions of HF and (NH4)2S208. 

The last upgrade concerned the temperature measure¬ 
ment. Since, the thermal expansion coefficient of ^®Si 
is about 2.6 x 10“® K“^, the measured volumes of the 
Si balls and unit cell must refer to the same tempera¬ 
ture to within a sub-millikelvin accuracy. Absolute tem¬ 
perature measurements are not necessary, but the tern- 
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FIG. 2. Phase of the x-ray fringes when the analyser is locked 
to an integer optical interference-order. Top: noisy supply 
voltage of the optical-interferometer actuators. Bottom: up¬ 
graded voltage supply. The integration time of each measure¬ 
ment is 30 s. The 192 pm period of the x-ray fringes has been 
used to show the analyser displacement (right scale). 



perature measuring-chains must be linked. We carried 
out more accurate and sensitive measurements of the Pt- 
thermometer resistance and linked our fixed point cells 
with those used to calibrate the measurements of the ^®Si 
ball temperature. 


B. Measurement procedure 


The measurement equation is 


a = 


VSd-. 


220 


y/SmX 

2n 


( 2 ) 


where ^220 is the spacing of the {220} planes, ac¬ 
counts for the different spacings of the {100} and {220} 
planes, and n is the number of x-ray fringes in a step of 
m optical fringes having period A/2. 

In practice, ^220 is determined by comparing the peri¬ 
ods of the x-ray and optical fringes. This is done by mea¬ 
suring the x-ray fringe fraction at the ends of increasing 
steps toA/ 2, where m = 1, 10, 100, 1000, and 3570. We 
start from A/(2d22o) = n/m Ri 1385.95 and measure the 
fringe fractions at the step ends with an accuracy suffi¬ 
cient for predicting the integer number of fringes in the 
next step. Consequently, A/(2d22o) is updated at each 
step. Eventually, measurements were carried out over 48 
subsequent steps of 3570A/2, 0.95 mm each for a total 
scanning length of 46 mm. 

The least-squares method is applied to reconstruct the 
x-ray fringes and to determine their phases at the ends 
of each stepi^ Typical input data are 300 photon-counts 
over 100 ms time windows spaced by 4 pm; a typical sam¬ 
ple contains six x-ray fringes, covers 1.2 nm, and lasts 30 
s. Each c ?220 measurement is the average of about nine 
values collected in measurement cycles where the anal¬ 
yser is repeatedly moved back and forth along the se¬ 
lected step. The visibility of the x-ray fringes approached 



FIG. 3. Lattice spacing values measured on March 14 along 
the line shown in Fig. [T] x rays enter the analyser from the 
obverse face. All the values are extrapolated to 20 °C, but 
not corrected for the systematic errors listed in Table H] The 
bars indicate the uncertainty of the linear interpolation of the 
8 detector-pixel values giving the value at 21 mm from the 
crystal base. The linear strain is due to the thermal gradient 
originated by the optical power injected into the crystal (from 
left side of the figure); the best fit line is also shown. The 
red dots indicate the outliers that were excluded from the 
subsequent analyses. 


50% with a mean brilliance of 500 counts mm“^. 
The crystal temperature is simultaneously measured with 
sub-millikelvin sensitivity and accuracy so that each d 220 
value is extrapolated on-line to 20 °C. 


C. Raw data 

Measurements were made over 48 subsequent analyser 
steps, 0.95 mm long. At each position, the lattice spac¬ 
ing was measured in the eight detector pixels and the 8 
results were processed to obtain, by linear regressions, 
the lattice spacing values in 48 points of the horizontal 
line that is the continuation of the laser beam, at 21 mm 
from the analyser base (Fig.[T|). A typical result is shown 
in Fig. El 

The hgure shows a gradient of the lattice spacing; we 
discovered that it is correlated to the temperature gra¬ 
dient caused by the power, about 0.75 mW, injected by 
the laser beam in the analyser. The reduced thermal 
conductivity of the residual gas in the vacuum chamber 
- because of the otherwise desirable low pressure - con¬ 
tributed to worsening the problem. The way we coped 
with this problem is described in Sec. IIV II 

Figure 2] shows that, after the thermal strain is re¬ 
moved, the residuals and the outliers are repeatable from 
one measurement to the next - also if carried out after 
one month. The head-on (obverse) and inverted (reverse) 
arrangements correspond to the analyser crystal mounted 
as it was in the boule and in a reversed arrangement; 
after the reversal, the x rays cross the crystal in the op¬ 
posite direction. The residuals and outliers repeatability, 
the scatter larger than the one expected by statistics, the 
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FIG. 4. Comparison of the variations of the measured lattice¬ 
spacing values along the line shown in Fig. [T] and after the 
thermal strain has been removed. Top: obverse face of the 
analyser, February 02 (blue) and March 14 (red). Bottom: 
reverse face of the analyser, May 15 (blue) and June 05 (red). 


different residual and outlier observed in the head-on and 
inverted surveys, and additional tests made by shifting 
the x-ray and optical baselines suggest that the outliers 
and residuals are caused by the analyser surfaces. 


Apart from the outliers, that were again pinpointed, 
the profiles shown in Fig. 2] are different from those given 
by Massa et al.^ Although the support-point variabil¬ 
ity could partially explain the difference and in the past 
measurements the resolution and repeatability were not 
as good as today (we only spotted a correlation between 
the profiles taken with the same analyser orientation and 
no correlation between those taken with opposite orienta¬ 
tions), we suspect that the difference is real and this sub¬ 
stantiates our surface-effects allegation. A collaboration 
with the Leibniz Institutes fiir Oberflachenmodifizierung 
in Liepzig is under way to develop machining technolo¬ 
gies based on plasma etching and ion beams to gain a 
better control of the geometrical, physical, and chemical 
properties of the crystal surfaces 


IV. ANALYSIS OF THE ERROR BUDGET 

A. Statistics 

Each measurement is the mean, after eliminating the 
outliers (red in Fig. [3]) and the thermal strain, of the sur¬ 
vey results. The uncertainty of the mean is dominated 
by the variations of the measured values, that are sup¬ 
posed to be caused by local effects of the analyser surface. 
Therefore, when calculating the uncertainty, we took the 
residual correlation into account. 


B. Laser beam wavelength 

The frequency of the Nd:YAG line is locked to the com¬ 
ponent R(56) of the 32-0 transition of the molecule. 
It was measured to better than a 10“^° relative uncer¬ 
tainty; therefore, it does not contribute to the measure¬ 
ment uncertainty. To eliminate the influence of the re¬ 
fractive index of air, the experiment is carried out in 
vacuo. With respect to our past measurements, the resid¬ 
ual pressure has been reduced by a factor of ten, to less 
than 0.01 Pa. Since the air refractivity at the atmo¬ 
spheric pressure is 2.9 x 10“^, assuming that the pressure 
is in the interval from zero to 0.01 Pa with a uniform 
probability, the relevant correction is 0.015(9) nm/m. 


C. Laser beam diffraction 

The period of the interference fringes is not equal to 
the plane-wave wavelength. In the case of the interference 
of two identical paraxial beams - whose angular spectra 
are strongly concentrated around a wave-vector having 
uj/c modulus, the difference between the period of the 
integrated interference pattern, Ag, and the plane-wave 
wavelength, A = 27rc/w, 

Ae - A _ (1 - xl/wl)Tv{V) ^ o? ^ 

where the optical-path difference of the interferometer 
arms is assumed to be much smaller than the Rayleigh 
length, xq is the offset between the beam axes measured 
at the beam waists, wq is the 1 /e^ spot radius at the 
beam waist, F is the central second-moment matrix of 
the angular power-spectrum of the beams, 2 q: is the mis¬ 
alignment between the beam axes, and 7 is the beam 
deviation from a normal incidence on the analyser. 

The Tr(r) term originates from diffraction and de¬ 
pends on the spread of the transverse impulse of the 
photons. It holds for any paraxial beam, no matter its 
profile is Gaussian or not^^; actually, it is obtained under 
a coaxial-beam assumption, i.e., when xq/wq = a = 0, 
and subsequently generalized to non-coaxial beams, but 
under a Gaussian-beam assumptions^ It must be noted 
that, in the case of Gaussian beams having cylindrical 
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Tr(r) / mrad^ 


FIG. 5. Radial profile of the focal plane image of the laser 
beam. The solid lines are minimum and maximum values 
of the radial profile of the best bivariate Gaussian function 
fitting the data. 


symmetry, Tr(r)/2 = 0 q/4, where 6*0 is the far-held di¬ 
vergence. We measured the angular power-spectrum of 
the beams emerging from the interferometer by using the 
Fourier transforming properties of a lens. Next, m is 
calculated from the central second-moment matrix of the 
focal plane image, which is recorded by a videocamera. 

Owing to the 8 bit resolution and dark-noise of the 
camera that we are presently using, a calculation of F 
based on a discrete approximation of the relevant inte¬ 
grals is unreliablei^ Therefore, it was estimated by ht- 
ting a bivariate Gaussian function to the focal-plane im¬ 
age; an example is shown in Fig. [S] The uncertainty 
associated to the Tr(r) estimate is small, typically, less 
than 1 %. 

To check the correction estimate, we examined the re¬ 
sults of a number of (^220 measurements carried out from 
2010 to 2014 with different beams. The results shown in 
Fig. [S] suggest that we overestimate the correction. Sub¬ 
sequent investigations did not shed light on this problem, 
but a study of the interference of wavefronts differently 
perturbed in the separate arms of the interferometer - 
where m does not hold exactly - seems to support an 
overestimation; more details will be given in a separate 
paper. Another hypothesis is a wrong estimate of the 
center of mass of the focal-plane image, which implies a 
correction always larger than true. In addition, since a 
single datum - corresponding to the largest beam diver¬ 
gence, see Fig. [5]- dictates the regression line, we may be 
mislead by a measurement error. Owing to the smallest 
beam divergence in the present set-up, the result we are 
reporting agrees with the value extrapolated to a zero 
correction from the data in Fig. [3] Therefore, we did 
not corrected the Tr(r)/2 = 3.87(58) nrad^ value; but, 
cautiously, increased its uncertainty to 15%. 

The interfering beams are kept parallel to within a 
2 a = 1 fivad maximum misalignment by levelling the 
phase in four quadrants of the interference pattern via the 
piezoelectric supports (pitch) and inertial drivers (yaw) 
of the interferometer base-plate. Consequently, the a^/2 


FIG. 6. Gorrected d 220 values vs. the applied correction for 
the diffraction of the laser beam. The measurements are made 
from 2010 to 2014 by using laser beams differently collimated. 
Bullets (red): He-Ne laser source @ 633 nm; squares (green): 
frequency doubled Nd:YAG laser source @ 532 nm. If the 
corrections are correctly estimated, the data are expected to 
lie on an horizontal line. 

term in m is irrelevant and was omitted. 


D. Laser beam alignment 

When assembling the apparatus, the laser beam devia¬ 
tion 7 from a normal incidence on the analyser was nulli¬ 
fied with the aid of an autocollimator looking at both the 
analyser and beam from the interferometer output-port. 
Next, as shown in Fig.[7l we carried out a number of d 22 o 
measurements while 7 was purposely changed along two 
orthogonal directions and its variations were recorded by 
an on-line telescope. The telescope is mounted, inside 
the vacuum chamber, on the same base plate of the x- 
ray/optical interferometer and picks up part of the out¬ 
put beam. After two parabola were fitted to the mea¬ 
sured ^220 values, the beam direction corresponding to 
the maxima - hence, to a supposed normal incidence 
- was identified to within ^rad uncertainty and main¬ 
tained in the telescope optics with an uncertainty of 20 
/rrad. Eventually, the laser beam was kept parallel to 
that direction and the telescope readings recorded for 
subsequent analyses. 

After we completed the measurements and removed 
the interferometer from the apparatus, we realised that 
the operation of the optical interferometer may be liable 
to a systematic error. This problem will be examined in 
a separate paper, but, since it relates to the assessment 
of the measurement uncertainty, we outline it shortly. 

In the case of a pointing error 7 , an analyser displace¬ 
ment s shears the interfering beams by 2sj. Hence, the 
measure beam goes through different parts of the optics 
crossed in its way to the detector. This shear changes the 
optical-path length by {2js)l3An, where An and P are 
the refractivity and the relevant component of the verti¬ 
cal angle of a wedge that, in a simplified model, substi- 
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FIG. 7. Measured d 220 values vs. the pitch angle of the laser 
beam; the angle origin is set in the maximum of the parabola 
(red line) that best fit the data. 


tutes for the optics in the way from the analyser to the 
detector. In addition, because of the wavefront curva¬ 
ture, the beam shear is sensed by the interferometer as a 
rotation equal to 



(4) 


where 1/i? is the wavefront curvature. Therefore, the 
measured ^220 value is 


dm = d: 


220 


1-V 
2 


f3An 


2R 


(5) 


where b is the Abbe’s offset between the laser and x-ray 
beam-centroids. According to ([5]), dm is maximum when 
7 = /3An -I- b/{2R), not when 7 = 0 as assumed in the 
alignment procedure. It must be noted that, when dm is 
maximum, the sensed rotation is not zero - as expected 
if 7 = 0 rad - but, 


^ sPAn 
^ R 


( 6 ) 


where, for the sake of simplicity, we assumed 6 = 0 mm. 

Since the analyser attitude is servoed so as to nullify 
signal of the angle interferometer - the differential phase 
between the quadrants of the interference pattern, the 
shear is counteracted by an analyser rotation. Eventu¬ 
ally, the pitch component of this rotation is disclosed by 
a (^220 gradient in the different detector pixels. The pitch 
explaining the gradient observed with a varying align¬ 
ment of the laser beam is shown in Fig. |8l When the 
beam is aligned in such a way that the ^220 value is max¬ 
imum, the pitch was equal to 0.4 nrad/mm in February 
and, after the analyser reversal and realignment, to 1.3 
nrad/mm in May. The yaw rotation might be similar, 
but, presently, we cannot detect it. 

Shear strains of the crystal lattice and surface effects 
(see Sec. IIII Cl) mimic the same gradients. Therefore, 
we cannot unambiguously explain a (6220 gradient by a 
parasitic pitch rotation. However, a number of additional 


FIG. 8 . Parasitic pitch of the analyser motion vs. the pitch 
angle of the laser beam. The analyser motion is servoed so 
as to nullify the differential signals detected by the optical in¬ 
terferometer. Measurements were done in February (squares, 
obverse face) and May (bullets, reverse face). The lines (solid 
and dashed, respectively) that best fit the data are also shown. 
The angle origin is set in the maximum of the parabola that 
best fit the measured ^220 values (see Fig. [3). 


tests excluded large strains. For instance, Fig. [5] shows 
the parasitic pitch of the analyser that explains the c ?220 
gradients observed in the March 14 and May 15 surveys, 
where the laser beam was aligned in such a way that 
the d 22 o measure is maximum. Since the only survey 
difference is the reversed analyser alignment, if we had 
observed a lattice strain, the two plots should be similar. 
Contrary, they are not; the mean pitch was 1.0 nrad/mm 
in March and —0.5 nrad/mm in May. 

According to (|1]), the radius of curvature of the inter¬ 
fering wavefronts, R = 25(1) m, can be estimated from 
the mean slope, s/R = 40(2) x 10“® mm“^, of the lines 
that best ht the data in Fig. [51 By explaining the (6220 
gradients in terms of parasitic pitch rotations and m, 
we estimate that the pitch component of /3 is 50 /irad 
from the March data and —25 firad from the May one. 
Since no change was made in the optical interferometer, 
these contradictory estimates exclude wedge angles much 
greater than, say, 50 ^rad. 

In order to estimate the needed correction and the as¬ 
sociated uncertainty, we used An Ri 1/2 and calculated 
the mean and variance of the c-y = (7 — P)^d 22 o/ 2 . cor¬ 
rection, where /? = 0(50) prad and 7 = 0.5/3 ± 20 /rrad 
are independently normally distributed. The final result 
is c-y = —0.11(48) nrad^. 


E. Laser beam walk 

The beam walk refers to the transverse motion of the 
interfering beams through the optical components. It 
originates from different effects causing the beams to 
move across imperfect surfaces or wedged optics. The 
effect of walks caused by a tilted incidence on the anal¬ 
yser was investigated in Sec. IIVDI In our previous set¬ 
up, the beam-splitter imperfection, combined with tilts of 


























7 




FIG. 9. Parasitic pitch of the analyser motion measured over 
0.95 steps in 48 subseqnent analyser positions. Surveys were 
carried out on February 02 (bottom, obverse face) and May 
15 (top, reverse face). The laser beam was pointed so as to 
ensnre that the measured ^220 value was maximum and the 
analyser was servoed so as to nullify the differential signals 
detected by the optical interferometer. 


the apparatus baseplate with respect to the laser beam, 
caused systematic differential variations of the optical 
paths through the interferometer that required ad hoc 
corrections. In the new apparatus, we made this prob¬ 
lem harmless by using a plate beam-splitter, having a 
parallelism error less than 10 ^rad, and by controlling 
electronically the baseplate level and tilt to within 25 
nm and 70 nm/m, over any analyser (short or long) dis¬ 
placement. The differential beam walk due to analyser 
parasitic rotations is irrelevant because rotations are less 
than 1 nrad/mm (see Sec. IIVDI) and the detector dis¬ 
tance is less than 0.5 m. 

The mechanical load driving the analyser carriage does 
the apparatus to sag and to yaw with respect to the laser 
beam. The relevant beam walks are quite large, up to 1 
fim and 5 p,rad, but, no correlation with the carriage dis¬ 
placement was observed. In addition, after any displace¬ 
ment, the link between the apparatus and the driving 
system is removed, thus allowing the same equilibrium 
position to be restored. It is difficult to estimate if there 
is a residual systematic differential-walk of the interfer¬ 
ing beams. If, over I mm analyser step, the systematic 
walk is in the [—0.1,0.!] mm interval, a 10% of what 
observed, and the differential wedge-angle between the 
end surfaces of the separate paths through the interfer¬ 
ometer is in the [—10,10] ^rad interval, the beam walk 
contribution to the uncertainty budget is 0.577 nm/m. 


F. Abbe’s error 


FIG. 10. Vertical variations of the x-ray fringe phase when the 
analyser is rotated - while keeping the displacement measured 
by the optical interferometer null - about the y (solid line) 
and z (dashed line) axes. 



FIG. 11. Vertical variations of the d 220 value measured over 
different 0.95 mm steps. The coordinate of the step center is 
also given. The solid lines are the linear regressions used to 
interpolate the measured values. 



The Abbe error refers to the difference, s • (b x fJ) = 
b • (J1 X s) = rj • (s X b), of the displacements sensed 
by the laser and x-ray interferometers, where b is the 
interferometer offset and Jl and s are the rotation and 
movement-direction of the analyser. 

As regards as H x s, it was zeroed to within 1 nrad 


FIG. 12. Solid line: variance of the residuals from the line 
best htting the d 220 values interpolated at different detector 
pixels. Dashed line: correlation between the same residuals 
and the residuals from the line best fitting the analyser pitch. 
The solid lines are expected variance and correlation. 



































(see Sec. lIV Dl) by servoing the motion so as the signals of 
the angle interferometer - the phase differences between 
the vertical and horizontal quadrants of the interference 
pattern - are null. 

As regards as s x b, the vertical offset was nullified by 
carrying out off-line measurements of the variations of 
the x-ray fringe phase in different detector pixels while 
the pitch component of ft is purposely changed while 
keeping the analyser displacement null. As shown in Fig. 
m we identihed the virtual pixel having a zero offset to 
within a 0.1 mm uncertainty. The horizontal offset was 
set to zero to within the same uncertainty by rotating the 
analyser about the vertical and by shifting horizontally 
the laser beam up to no phase variation is detected, as 
shown in Fig. [TUI 

The angle interferometer and attitude control display 
imperfections; we took advantage of the resulting pitch 
noise - which is shown in Fig. 0 - to check the vertical 
offset by data analysis. The eight detector-pixels have 
a linearly increasing offset and, as shown in Fig. [TTJ the 
linear regressions of the ^220 values intersect, ideally, in 
the pixel having a null offset. Since the ^220 value in this 
pixel is insensitive to the pitch noise, the best way to 
find it is to look at the minimum variance of the resid¬ 
uals from the best-fit lines of the values interpolated in 
each detector pixel or, which is the same, at the zero- 
correlation between the same residuals and the residuals 
from the best-fit line of the pitch noise. Examples of 
the residual variance and correlation are shown in Fig. 
[T2I In half of the surveys the null offset locates in a 
pixel that differs from where the phase variation of the 
x-ray fringe was found - in previous off-line experiments 
- insensitive to the analyser pitch rotation. After the 
projection on the analyser, the maximum differences are 
0.5 mm. Since the x-ray/optical interferometer rests on 
an double anti-vibration system (a table supported by 
air springs, in turn, mounted on a giant pendulum), we 
explained these shifts by instabilities of the relative lev¬ 
elling between the x-ray source and interferometer which 
affect the pixel looking at the analyser point having zero 
vertical offset. 

We trusted the zero-offset pixel identihed by the data 
analysis. Since we carried out this analysis after the in¬ 
terferometer was removed from the apparatus and, con¬ 
sequently, we did not investigated experimentally the 
problem, we assumed a null-pixel error uniform in the 
[—0.5, -1-0.5] mm interval with an uncertainty of 0.29 mm. 
By combining this uncertainty with the uncertainty of 
the zeroing of the horizontal offset component, 0.1 mm, 
and parasitic rotation, 2 nrad, we estimated that the 
Abbe error was nullihed to within a total uncertainty 
of 6.11 X 10-i%20- 


G. Movement direction 

The analyser moves orthogonally to the front mirror 
of a trihedron; straightness errors are nullihed to within 



FIG. 13. Relative directions of the reciprocal vector h of 
the diffracting planes, unit normal to the analyser n (obverse 
mirror), and displacement direction s. The dates indicate 
the s checks; after 2014/4/14 the analyser was removed and 
realigned in reverse. The circles indicate the measurement 
uncertainties. The dashed line is the locus of the s direction 
where the projection error is null. 


nanometers by servoing the motion with the signals of 
capacitive transducers that sense the transverse displace¬ 
ments of the trihedron top- and side-face. 

The misalignment between the optical and x-ray inter¬ 
ferometers causes them to measure different components 
of the displacement. With s indicating the displacement, 
the x-ray interferometer senses s • h, where h is the unit 
normal to the diffracting planes, whereas the optical in¬ 
terferometer senses s • n, where n is the unit normal to 
the analyser, see Sec. IlVDI The difference, s • (n — h), is 
null when the displacement is orthogonal to (n —h), that 
is, when it bisects the angle formed by h and h. The 
lattice constant is linked to the measured mA/(2n) ratio 
by 


<^220 — 


mX 

2n 


1 -I- s • 


(n 



(7) 


The two (obverse and reverse) analyser-mirrors are pol¬ 
ished parallel to the {220} planes. The residual misalign¬ 
ments, jn — h| = 13.7(1.2) /rrad and jn —h| = 10.8(1.2) 
/rrad for the obverse and reverse mirrors, respectively, 
were estimated by a least-squares adjustment of the mis¬ 
alignment between x-ray and light reflections on the mir¬ 
rors and lattice planes^ the phase shift of the x-ray 
fringes when the analyser motion lies in the mirror planes, 
and the measured angle between the two mirrors. In or¬ 
der to calculate the relevant correction and the associated 
uncertainty, the (s —n) angle - that is, the angle between 
the trihedron and analyser - was periodically measured 
to within a 10 /irad uncertainty; examples of the mea¬ 
surement results are shown in Fig. 1131 
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H. Analyser temperature 

The analyser temperature is measured by a capsule 
standard Pt resistance thermometer inserted into a well 
in a copper block in thermal contact with the crystal. 
Resistance measurement were carried out by a FLUKE 
1595A multimeter; to minimize the self-heating, the mea¬ 
surement current was 0.3 mA. Each temperature datum 
is the average of 15 measurements pairs, carried out with 
both positive and negative currents and integrated over 
30 s. The 1595A linearity was checked by a resistor net¬ 
work made in such a way that the voltages across any 
number of resistors in a resistor series are read to get 
four-terminal values interrelated by the formula for the 
series connection^ The test showed that linearity is bet¬ 
ter than 10 /iU - corresponding to 25 /rK - for resistance 
measurements from 90 to 120 U. 

The lattice constant measurements were carried out 
from February to June 2014; on April 08, the thermome¬ 
ter was calibrated in situ ~ that is, by moving the ther¬ 
mometer from the apparatus to the fixed-point cells with¬ 
out changing the measuring chain and cables. We ex¬ 
trapolated the resistance readings to a zero current and 
corrected the cell temperatures for the immersion depth 
and hydrostatic pressure. 

The temperature measurements require sub-mK accu¬ 
racies and any difference between the temperature scales 
used to extrapolate the molar volume and lattice con¬ 
stant to 20 °C must be excluded or identified. Conse¬ 
quently, on April, our fixed-point cells were compared 
with those of the Physikalish Technische Bundesanstalt. 
After the corrections for the immersion depth and hydro¬ 
static pressure were taken into account, the differences of 
the resistance readings were 

Rtpw(PTB) - Rtpw(INRIM) = 5(3) 

RGa(PTB) - i?Ga(INRlM) = -7(4) 

Unfortunately, it was not possible to investigate the non¬ 
uniqueness associated with the readings of the two ther¬ 
mometers at 20 °C; it was cautiously set to 0.1 mK.— 

Taken note of these differences, the uncertainties of the 
triple point of water and melting point of Ga realisations 
are irrelevant; the repeatability of the cells is 50 /rK. Ow¬ 
ing to the huge data averaging, the noise of the resistance 
measurements is irrelevant; the stability of the reference 
100 U resistor over the two months before and after the 
calibration is 33 /rfl, the linearity of the measurement of 
the 107 O thermometer-resistance is better than 10 /rO, 
the measurement non-uniqueness is 0.1 mK. All together, 
the uncertainty of the temperature measurements, esti¬ 
mated by Monte Carlo simulation, is 0.17 mK. 

Each 6,220 measurement was extrapolated to 20 °C ac¬ 
cording to2^ 

d22o(7o) = d22o{T) [l -I- ai(ro — T) + a2{To — T)^] , 

( 8 ) 

where Tq = 20 °C, ai = 2.5530(12) x 10"® K'^ and 
a 2 = 4.32(37) x 10“® K~^. All measurements were car¬ 
ried out in the temperature range from 19.9 °C to 20.3 



FIG. 14. Comparison of the finite-element calculation of the 
minimum (the thermal flux is grounded through the ther¬ 
mometer copper block) and maximum (the thermal flux is 
grounded through the crystal support points) thermal strains 
(with respect to the the measured crystal temperature) due 
to the optical power, 0.75 mW, injected into the analyser by 
the laser beam (solid lines, red) and the least-squares adjust¬ 
ment (filled area) of the ^220 gradients and variations at the 
crystal ends (dots). 


°C; therefore, the average extrapolation uncertainty is 
0.24 X 10-9^220- 

The thermometer self-heating was identified by repeat¬ 
ing 6,220 measurements with varying currents; the rel¬ 
evant correction for the 0.3 mA current is 0.33(3) x 
10 “9(i220i the measured value being smaller than the true 
one. 

The calibration history, dating back to December 2007, 
shows a linear drift of 14(5) /rU/month or 0.035(13) 
mK/month that was taken into account to extrapolate 
the calibration to the actual measurement date. 

Eventually, the total uncertainty of the lattice constant 
extrapolation to 20 °C is 0.497 x 10“9d22o- 


I. Thermal strain 


A linear approximation of thermal strain due to the op¬ 
tical power injected into the analyser by the laser beam. 


^6220 

6220 


= a — b(x — xo), 


(9) 


where a = 0.948(203) x lO"®, b = 0.036(9) x lO'^ mm"!, 
and xq = 22.8 mm is the center of the survey, was found 
by a least-squares adjustment of the 6220 gradients and 
the results of repeated 6220 measurements carried out 
with varying optical powers. The comparison of (jH)) with 
the numerical calculations of the thermal strain is shown 
in Fig. [TT] To correct for the thermal strain, we trusted 
the value given by (0, but increased its uncertainty to 
the one half of the 1.28 x 10“® gap between the minimum 
and maximum strain predicted by the numerical calcula¬ 
tion. Therefore, the interpolated value at 22.8 mm was 
reduced by 0.948(640) x 10“9d22o- More details about the 
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FIG. 15. Residual mean self-weigh strain of the analyser cal¬ 
culated at 21 mm from the base as a function of the distance 
between the support points. The filled curve is the distance 
probability-distribution, given three contact points uniformly 
distributed in the (5 x 5) mm^ support areas. 

numerical and experimental investigations of the analyser 
response to the thermal load will be given in a separate 
paper. 


J. Self-weight deformation 

To average the ^220 measurements over the largest 
crystal part, a long analyser has been used. The sim¬ 
ulation of the gravitational bending allowed the analyser 
to be optimally designed, the residual lattice strain to 
be predicted, and the contribution of the self-weight de¬ 
formation to the uncertainty budget estimatedi ^°’^^ The 
simulation purposes were to find the maximum height of 
the analyser lamella consistent with a non-strained lat¬ 
tice and the support points minimizing bending or sag¬ 
ging. In order to estimate the necessary correction and 
its uncertainty, the residual strain at the 20.5 mm height 
was recalculated for points randomly located inside the 
(5 X 5) mm^ support areas. The results show that the 
mean strain, which is shown in Fig. 1151 depends only 
on the distance between the support points. Since the 
simulation indicates that the everted-inverted transition 
occurs with a slightly expanded lattice, we reduced the 
measured values by 0.543(377) x 10“®(i22o “ which relates 
with the distance distribution shown in Fig. 1151 


FIG. 16. Measured (£220 values. Each measurement is the 
mean, after eliminating the outliers and correcting for the 
thermal strain, of a survey of a 48 mm long crystal part. 

respect to the diffracting planes - and of local surface 
strain, if any, were washed out by the survey averaging. 

The linear gradient of the mean analyser thickness over 
the 48 mm measurement distance is less than 10 /im and 
contribute by less than 0.4 x 10“® to the ^220 measure¬ 
ment. This error is nullified by repeating the measure¬ 
ment after a 180° rotation of the analyser and by aver¬ 
aging the results. If the crystal displacement does not lie 
in the mean surface of the analyser, the interferometer 
defocuses. The out-of-plane angle is less than 2 /im/cm, 
to which a zero-mean uniform error having 0.23 nm/m 
standard deviation will correspond. 

A stress exists in the crystal surfaces even if the bulk 
material is stress-free. This problem was investigated by 
Quagliotti et. by using an elastic-film model to pro¬ 
vide a surface load in a finite element analysis. The study 
showed that, if the tensile stress is 1 N/m, the measured 
lattice spacing is 6 x 10“®d22o smaller than the value 
in an unstrained crystal. Literature values of the (001) 
surface-stress obtained from ab initio and molecular dy¬ 
namics calculations are given by Quagliotti et. at^ the 
stress of the (110) surface is expected to be 60% smaller. 
Owing to the value and sign scatters of the literature 
data, we do not propose a correction and associate to a 
null stress an uncertainty of 0.1 N/m. Therefore, the rel¬ 
evant contribution to the lattice constant uncertainty is 
0.6 X 10“®d22o- Further experimental investigations and 
atomistic calculations are under way to confirm that sur¬ 
face stress effects are irrelevant or to quantify and correct 
for themi^ 


K. Aberrations of the x-ray interferometer 

Geometric aberrations contribute to the phase of the 
x-ray fringes by less than 0.01c?22o per 1 ^m changes of 
the analyser thickness or focusingj ^^i^^ The root-mean- 
square roughness of the analyser surfaces is less than 1 
/im, with its main components in the neighbour of the 
0.1 mm wavelength. The effect of the surface roughness 
- included the large local variations of the angle with 


V. MEASUREMENT RESULTS 

Three ^220 surveys were made on February 12 and 
March 10 and 14 with the analyser in the head-on ar¬ 
rangement; four were made on May 14 and 15 and June 
05 and 12 with the analyser in the back orientation. Ex¬ 
amples are given in Figs. [3] and |H Next, after eliminating 
the outliers and correcting for the thermal strain, each 
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this work 


ref. [8] 


711.5 712.0 712.5 713.0 713.5 

(^220 - 192014000 am) / am 

FIG. 17. Comparison of the d 220 values GQl) and m- 

^220 profile was averaged to obtain the mean lattice spac¬ 
ing. 

The results are shown in Fig. [1^1 an example of the 
error budget is given in Table HI With respect to our pre¬ 
vious measurement, in addition to the reduction of the 
total uncertainty, the Table H] shows a significant redis¬ 
tribution of the uncertainty contributions. The values 
belonging to each obverse/reverse faces set are signifi¬ 
cantly correlated; this explains the repeatability, which 
is much better than the uncertainty. 

The analyser reversal required a full realignment of the 
two interferometers. Therefore, only the wavelength and 
temperature uncertainties, laser-beam diffraction, self- 
weigh deformation, and aberrations of the x-ray interfer¬ 
ometer combine in the same way; the remaining contri¬ 
butions to the total uncertainty are largely independent. 

Figure [TBl shows a difference between the values mea¬ 
sured with the analyser mounted in the head-on and in¬ 
verted arrangements. We are not yet able to say if this 
difference is real, e.g., due to a different physical and/or 
chemical structure of the analyser surfaces, or if it in¬ 
dicates that the control of the systematic errors is less 
good than what we estimated. The first option is sup¬ 
ported by the repeatable observation of different obverse- 
and reverse-profile, which was statistically anticipated by 
Massa et al.^ The second is supported by the fact that a 
difference between the obverse and reverse mean-values 
of d 22 o was not reported. 

The final measured value, 

d22o(2014) = 192014711.98(34) am, (10) 

at 20 °C and 0 Pa, is the mean of the data in Fig. [1^1 
The relative uncertainty is 1.75 nm/m. In the average, 
we did not take the data uncertainty and correlation into 
account, but, to avoid that the different number of ob¬ 
verse/reverse surveys biases the result, firstly, we aver¬ 
aged the obverse and reverse data and, subsequently, av¬ 
eraged the two results. To be conservative, we associated 
to the mean the worst uncertainty of the input data. 

As shown in Fig. [T71 (ITU)) is slightly smaller than 

( 11 ) 


TABLE I. Relative correction and uncertainty, in parts per 
10®, of the 2014/02/12 d 220 value. 


Contribution 

Correction 

Uncertainty 

data averaging 

0.000 

0.722 

wavelength 

-0.015 

0.009 

laser beam diffraction 

3.978 

0.597 

laser beam alignment 

-0.110 

0.480 

beam walks 

0.000 

0.577 

Abbe’s errors 

0.000 

0.611 

movement direction 

0.699 

0.214 

temperature 

-0.500 

0.497 

thermal strain 

-0.948 

0.641 

self-weigh 

-0.543 

0.377 

aberratioirs 

0.000 

0.642 

total 

2.56 

1.75 


given by Massa et a74 A reasons might be a positive bias 
of the correction for diffraction applied to (ED, as shown 
in Fig.ini In addition, in ED, the pointing stability of the 
laser beam was not monitored on-line and, in retrospect, 
we might have overestimated the pointing error. Since a 
pointing error causes a measurement-underestimate, 
we consistently - but, perhaps incorrectly - applied a 
relatively large positive correction. Eventually, the re¬ 
analysis of the effect of a non-orthogonal incidence of the 
laser beam on the analyser mirror in Sec. lIVDl shows that 
wedged optics in the beam path combine with a wrong 
pointing to originate a positive error. Therefore, contrary 
to what we did in the past, the correction for the laser 
beam alignment in Table H] is negative. 

To evaluate the correlation between ED and ED, the 
following considerations are in hands. The optical in¬ 
terferometer used for the two measurements were com¬ 
pletely different as regards the wavelength, laser source 
beam geometry and power, and optical components; the 
combined x-ray and optical system was also completely 
knocked down and realigned. Furthermore, the temper¬ 
ature measurements relied on a different, of higher qual¬ 
ity, apparatus and were calibrated anew. Eventually, all 
the contributions to the error budget were re-examined; 
some - the thermal strain, self-weigh deformation, sur¬ 
face stress, and walks of the laser beam -- were reconsid¬ 
ered or calculated anew. 


VI. CONCLUSIONS 

We re-checked the measurement uncertainty of the 
lattice parameter of the ^®Si crystal used to determine 
the Avogadro constant. Additional measurements and 
stress-tests were carried out by using of an upgraded mea¬ 
surement apparatus. No error was identified; this work 
conhrms the value given by Massa et al^ and yields an 
additional results having a reduced uncertainty. 


d22o(2011) = 192014712.67(67) am 
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